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^\ • We discuss a dynamical mechanism of pion production from disoriented chiral condensates (DCC) . 

0^ ', It leads to an explosive production of pions via the parametric resonance mechanism, which is similar 

to the reheating mechanism in inflationary cosmology. Classically it is related with the instability in 
^ ' the solutions of the Mathieu equation and we explore the quantum aspects of the mechanism. We 

O , show that nonlinearities and back reactions can be ignorable for sufflciently long time under the small 

'"T' ■ amplitude approximations of background a oscillations, which may be appropriate for the late stage 

of nonequilibrium phase transition. It allows us to obtain an explicit quantum state of the produced 

0^ , pions and a, the squeezed state of BCS type. Single particle distributions and two-pion correlation 

functions are computed within these approximations. The results obtained illuminate the charac- 
teristic features of multi-pion states produced through the parametric ampliflcation mechanism. In 
particular, two-pion correlations of various charge combinations contain back-to-back correlations 
which cannot be masked by the identical particle interference effect. We suggest that the parametric 

(^— N [ resonance mechanism might be a cause of the long lasting amplification of low momentum modes 

pf-\ . in linear sigma model simulations. 
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I. INTRODUCTION 

Oh! 



rr^' We now believe in quantum chromodynamics (QCD) as the theory of strong interaction. Yet, the Centauro and 

(-H ' anti-Centauro events that were found in cosmic ray experiments ||l| indicate some puzzling features. The events are 

• • characterized by large fluctuations in the ratio of number of neutral pions to that of charged pions. It appeared to be 

,J^ unlikely that such feature can be realized inside the conventional picture of hadronic interactions with "soft" QCD 

K^ , interaction, which implies more or less an independent emission of particles. 

^ ■ However, several authors formulated a scenario which explains such features of the events based on the idea of 

5t ! formation of chirally misaligned domains, the disoriented chiral condensate (DCC) 0,0]. It requires that a "hot" 
matter forms during hadronic collisions in which the chiral symmetry is restored, and then it rapidly cools down so 
that chiral orientation of the pion fields can align to a random direction different from that of the vacuum. If the low 
momentum mode of the pion fields are enhanced large domains of definite (but random in direction) isospin would 
form. Then, a large isospin, and therefore, charge fluctuations naturally results. 

The amplification of low momentum pion modes was demonstrated by Rajagopal and Wilczek B by doing numerical 
simulation of the linear sigma model using the quench initial condition. A more explicit confirmation of the formation 
of large domains was provided by Asakawa, Huang and Wang H who computed pion-field correlation functions by 
using an elaborated simulation code which takes into account expansion along the longitudinal as well as the transverse 
direction. 

One should note, however, that the large charge-neutral fiuctuation does not necessarily imply a random isotropic 
in isospin space) rolling down from the top of the Mexican hat potential. As was shown by Greiner, Gong and Miiller 
any isosinglet multipion states approximately have the well known neutral fraction distribution P{f) 



m-^. (1) 



where 



N^ 

We note that this distribution was in fact derived by only requiring the isospin singlet phase space [Q. Therefore, the 
large charge-neutral fluctuation may merely imply isospin singlet DCC state, not a random isotropic rolling down. 

The interpretation of the results of the simulation by Rajagopal and Wilczek B also requires reexamination. They 
interpreted the origin of the amplification of the low momentum pion mode as due to instability of the Nambu- 
Goldstone modes of chirally symmetric field configurations at around the top of the Mexican hat. However, the 
enhancement of the low momentum pion modes actually take place with much longer time scale than that of rolling 
down motion, ~ l/m„, as one can clearly see in Fig. 1 in their paper. Then, we need to identify certain mechanism 
which is responsible for the enhancement. The most likely candidate, to our opinion, is the parametric resonance 
mechanism [g|,y] . The mechanism, possibly being tied up with a new instability that arises within the approximation 
scheme we will take, seems to be capable of explaining qualitative features of long lasting amplification of low mo- 
mentum modes. (See sections 3 and 6.) Some other aspects of formation mechanism of DCC were explored in Refs. 

In a previous paper |15| , we briefiy explored the parametric resonance mechanism emphasizing its characteristic 
features in the two pion correlations, and its possible relevance as a hunting tool of DCC. In this paper we examine 
the mechanism in detail, hoping that it will shed some light to the origin of enhancement of the low momentum 
pion modes. We focus on the pion production at around the bottom of the Mexican hat potential. It may be an 
appropriate setting at least for the late stage of non-equilibrium phase transition which would result in the formation 
of coherent DCC domains. 

Furthermore, we concentrate on the case that the sigma model fields oscillate along the a direction. In fact, it is 
shown in the numerical simulation g that sigma model field configurations rapidly point to the a direction in a time 
scale of the order of ~ I/ttIt^. Under the background oscillation of the a field the quantum fiuctuations of the pion 
and the a fields can be excited. If the frequencies of background and quantum fluctuations match with each other 
the fluctuations are parametrically amplified, leading to an explosive pion production. 

We emphasize that the mechanism of pion production from DCC which we explore in this paper differs from that 
is usually adopted in the literature &-H] as we mentioned earlier. In the conventional picture of pion emission from 
DCC, one assumes that the sigma model fields roll down along the direction (in most cases) different from the a 
direction. Then, the field configurations relax toward the a direction, the orientation of the QCD vacuum, which 
entails the coherent pion emission through the relaxation process. In contrast, pion production takes place even when 
the sigma model fields rolls down along the a direction in our parametric resonance mechanism. To indicate this 
point, we propose in the next section a set of approximations that leads to a simplified but a concrete model field 
theory that allows us to construct explicitly the produced multi-pion quantum state. 

We should remark that we do not claim that the parametric resonance is the whole story. Formation and decay 
of DCC is a complicated nonlinear phenomenon. The instability associated with the Nambu-Goldstone modes is 
indispensable ingredient for triggering the growth of the low momentum pion modes. But, we also stress that it is 
not enough either. In Sec. VI we will calculate the power spectrum of pion and sigma fields in our framework and 
attempt a qualitative comparison with the result obtained in Rajagopal-Wilczek's simulation. It will give us a feeling 
on the question of to what extent the amplification of power is due to the parametric resonance mechanism. 

Even though the parametric resonance mechanism is the cause of the long las ting amplification of low momentum 
modes we still lack understanding how the initial instability is mediated to it p4|| . Leaving understanding of this 
last point to the future investigations our emphasis in this paper is that if the parametric resonance mechanism is 
operative at the late stage of the non-equilibrium phase transition it may give clear signatures which should be useful 
for its experimental hunting. 

In Sec. II, we introduce the linear sigma model, formulate the system with parametric resonance by introducing 
successive approximations in various stages, and make clear that to what extent the approximations are valid. In Sec. 
Ill, we quantize the system in a conventional way of quantizing a system of quantum fields, and clarify the relationship 
with the other formulation. In Sec. IV, we calculate the single particle distributions of pions and a to demonstrate a 
characteristic feature of the parametric resonance mechanism. In Sec. V we discuss the two pion correlations which 
display another characteristic feature of the mechanism. We argue that they can be used as powerful experimental 
hunting tool of DCC. In Sec. VI we compute the power spectrum of pion and sigma fields and examine the qualitative 
features of the amplification of power in our framework. In the last section we will summarize our investigation, 
give concluding remarks, and discuss the limitations of our framework. In Appendix A, we summarize the feature of 
multiparticle states implied by the single mode squeezed state. In Appendix B, we derive sum rule obeyed by isospin 
invariance and the iso-singlet nature of the multiparticle state within the framework of single mode treatment. 



II. THE MODEL AND THE APPROXIMATIONS 

The Lagrangian of the linear sigma model is given by 

C^]^d^,d)ad^'cj)a-^{c^a<t>a-vlf+ha. (3) 

where (j)a = (ct, tt). The typical values of the parameters which are relevant for phenomenology of QCD are: 

A = 20, vo = 90MeV, m^ ^ J — ^ 140MeV. (4) 

V vo 

These values of the parameters will be used in the numerical analysis in later sections. 

As discussed in the previous section we deal with the pion production that takes place when the sigma model fields 
are at around the bottom of the Mexican hat potential, which may be realized in the time scale of '^ ^Ifn-^ after 
rolling down. We assume that at this stage the sigma model fields oscillate along the a direction. For clarity we 
construct a further simplified model by doing successive approximations in the linear sigma model that allow us to 
explicitly construct the produced multi-pion quantum state. 

We expand the sigma model fields at around the minimum of the Mexican hat potential along the a direction. 



X 



The linear sigma model Lagrangian, then, takes the form 



(5) 



\{d.xf~\{iv^-vl)x^-\vx'-\x'^\^ 



1 2 ^2 \ -<1 1 \ 2 ^2 -'^ \ ^4 

'2^1^'^ -AwxTT - -Ax TT - -Att . 



(6) 



Following Mrowczyhski and Miiller |9|], we make a further approximation of decomposing a model fields into the 
classical time-dependent background fields and the quantum fluctuations around them 

7r(a;,i) = 7ro(t) +^(a?,t). 

We have assumed the spatial homogeneity of the background fields. It may be a reasonable assumption if we work 
inside a single domain because the key feature of DCC is its coherence over a domain. Of course, it is a simplifying 
assumption in real physical situation, but is an inevitable one in order to make the treatment technically manageable. 
Substituting (0) into Lagrangian (jg), we have 

I^\xq + ?7, 7?o + C] = -Cixo, 7?o] 



(8) 



where rricr — •\/A(3f^ — Vq) and v is the minimum value of the a direction determined by the equation X{a'^ — VQ)a—h = 
0. For small h, v is given approximately by w = wq + 2\y2 ■ Due to the background field oscillation, a time-dependent 
term E(i) arises in the mass terms of ^ and rj fluctuations and is given by S = Axo(xo + 2w). The terms linearly 
proportional to fluctuations, of course, vanish due to the equations of motion of the background fields: 

Xo + Axg + 3Awxo + "I'xo + A(xo + w)^o = 0, (9) 



2^ 



TTo + m^TTQ + 2Awxo7ro + ^X^T^a + Attq = 



(10) 



We work with the Ansatz that the background fields osciUate along the a direction, and set ttq = in the subsequent 
analyses. We further restrict ourselves into the regime, | — | <C 1, in which the nonlinear terms in xo's equation of 
motion are negligible. Then, the solution to the xq's equation of motion takes the harmonic form 



Xo{i) == Xocos(m<^i + </?), 



(11) 



where (p is an initial phase which we set to vanish for simplicity. In the numerical analysis performed in Sec. IV, we 
will use | — I = 0.05 (x 1/2, 1, 2) for illustrative purpose. 

We should note that the limitation of small amplitude oscillation that we have imposed, j-^l ^ 1, restricts the 
validity of our treatment only to a qualitative level. Nevertheless, we believe that it is a meaningful starting point, 
and in fact it makes theoretical analysis transparent as we will see below. 

We differ from Mrowczyhski and Miiller [p| by treating 77 and ^ quantum mechanically. The equations of motion 
obeyed by fluctuations 77 and S, are given as 

d'^ri + mlr] + 3(2Awxo + >^Xo)'n + 2A(xo + v)v'^ + >^V^ = 0, 

(12) 
52^ + m^e + (2A«xo + Ax§)C + A^^ = 0. 

They are greatly simplified by the restriction of small oscillation amplitudes of background fields. 

Let us first focus on the pion fiuctuations. The Xo term is ignorable compared to vxo term because of the restriction 
of 1^1 ^ 1. The vxo term is also small compared with the mass term, but we shall keep it otherwise we lose the 

parametric resonance. The cubic term of f is negligible until the time scale 



4Awxo 



■In 



(^)'-2At;xo^ 



lOA 



(13) 



where the dimensionless time z is measured by m^ as 2z = rria-t + ip + TT. For instance, it can be numerically estimated 
as z- 67 for |^| =0.05. 

We have ignored the quantum back reactions to the pionic and sigma fluctuations due to particle production. It 
could be taken into account, e.g., by the Hartree-Fock approximation. Instead of going through this treatment, we 
estimate the time scale until which it can be negligible. It is given roughly as 



■In 



(. 



mlxa 



4Awxo \10A(xo 



(14) 



and the same numerical examination indicates that it is at 2; '--^ 32. 

We emphasize, therefore, that under the approximation of small amplitude oscillation of background fields, ignoring 
nonlinear terms in the field fiuctuations is a good approximation in a fairly long time even for such strong coupling 
as A = 20. 

We end up with the equation of motion of sigma and pionic fluctuations: 



-p2+ Aa- 2q„ cos(2z) 



4'^(^) = o, 



(15) 



r ^2 



rfz2 



where 



and 



A-^ — 2q-„ cos(2z) 



efe(^)=0, 



_ 4(fc2+m2) _ 12Ai;xo 



mt, 



mt, 



^ 4(fc2 + ml) _ AXvxQ 



(16) 



(17) 



(18) 



rjf^ and ^j^ denote the Fourier components of the fields r]x and ^{x), respectively. The equations (nS) and ( pq ) are 
known as the Mathieu equation and are known to admit unstable solutions for a wide range of parameters of A's 
and q's. See, e.g., Ref. p6[ . Such unstable solution may be interpreted as explosive particle production under the 
background of oscillating xo fields. 



III. QUANTUM EVOLUTION OF THE SYSTEM AND THE SQUEEZED STATE 



We discuss in this section the quantum evolution of the time-dependent states of 7? and a fluctuations. Thanks to 
the restrictions and approximations introduced in the previous section, the structure of quantum states formed by ry 
and ^ quanta can be analyzed analytically. 

We start by writing the Hamiltonian in a form analogous to the harmonic oscillators but with time-dependent 
frequencies: 



B 



J 



d^k 



\pkP-k+\m^fQkQ-k 



E/^^^[^^^^!fc + ^"feW^QiQifc 



(19) 



where 



ni{tf = k'+ml + &\vxo{t), 



The variable Q^ and Q-'i in (|9|) are defined by 



ml + 2Xvxo{t). 



(20) 



r]{x,t) 



(fk 

'{2^ 



ikx 



e^"-^QfcW, e(a=,t) = 



(fik 

(2^ 



e'fc-^gyo, 



(21) 



and Pf^ and P^ are, as usual, the conjugate momenta of Q^ and Q^., respectively. The index j runs over 1 to 3, and 
we take the adjoint representation for the iso-triplet pion fields. For its concrete form, we refer Appendix B. One can 
see that the equations of motion (15) and ( |l6|) are derived from this Hamiltonian. 

Within the small background oscillation we can take that (— ) < (2iiL)2^ This restriction guarantees that the 
frequency f^r._ is real, which is necessary in our present treatment which ignore quantum back reactions. If the 
background oscillation amplitude is larger than the critical value the system have an instability which we referred to 
as a new instability in Introduction. It definitely arises within our approximation and is indicative of a remnant of 
the Rajagopal-Wilczek instability because it arises only for negative xo (which means toward central maximum of the 
Mexican hat). Nevertheless, the new instability has a characteristic feature that is quite different from that of the 
Rajagopal-Wilczek's. The former can continue for long time while the latter lasts only a time scale of l/rria. Notice 
that both of the instabilities that we are talking about are the local instabilities, not the real global ones which totally 
destabilize the whole system. 

The nature of the instability, however, is not transparent to us. When one works with certain truncation in highly 
nonlinear systems one easily generate instabilities which may or may not be possessed by the original system. There 
are instabilities that one can easily understand their origin from simple consideration. For example, pion fields can 
have local instability by generating the winding motion encircling around the bottom of the wine bottle when the 
field energy of background a field oscillation with amplitude xo exceeds the symmetry breaking energy, that is, 
i?Ti^Xo > m^v^ . It leads to the condition for onset of the instability ^ > ^ 

But the instability we just encountered above occurs at much smaller amplitude, (— ) > (I2^)2_ r^^ ^^j. knowledge 
there is no intuitive way of understanding the instability. It can be the artifact of our approximation and truncation 
scheme and, if this is the case, may be cured by taking nonlinearities and back reactions into account. 

In the rest of the paper we will restrict ourselves into the small amplitude regime of the background a oscillation 
so that we are free from the new instability. In Sec. VII we will briefiy discuss its role together with the parametric 
resonance as a possible candidate mechanism for understanding the long-lasting amplification of powers in sigma 
model simulations. 

The quantum theory of harmonic oscillators with time-dependent frequencies were developed since sometime ago . 
In a previous paper jl|] we followed the formalism developed by Shtanov, Traschen and Brandenberger (STB) p7[ , 
but we shall present in this paper an equivalent but more conventional formalism of quantizing the same system. 
From the form of the Hamiltonian (HOh , we can discuss sigma and pion sector collectively; we treat below the sigma 
sector but with suppressing the superscript a, and it can be regarded as that of pion sector in which the superscript 
TT as well as iso-spin index are suppressed. 

We decompose the dynamical variables in momentum space in terms of solutions of the equation of motions 



where Qi and Q2 are the solutions to the equation of motion of Q derived from the Hamiltonian, 

gfc(t) + fi|(t)Qfc(i)-o. 



(22) 



(23) 



Notice that Q2 — Ql- a„i, and a , are the annihilation and the creation operators, respectively. They are time- 
independent, i.e., they obey the Heisenberg equation of motion as 



da 



ok _ ^%k 



+ ^[H,%k\ =0- 



dt dt 

The consistency condition between the commutation relations, i.e., 



J 



1' [Qfc, Pfc/] - i5(fc - fc') 



«0fc' %fc 

is achieved by the Wronskian condition 

Q1Q2 - Q1Q2 = i- 
Substituting ( p^ into (pj|), we obtain the following form 



H, 



k 



d^k 



1 



-.t J 



^k%k%-k + ^>ok%~k + ^^Kk%k + %k%k^ 



where the coefficients A^^ and B^^ are defined by 



Aj^ = Ql + nlQl, B,^^QiQ2 + nlQiQ2. 



(24) 



(25) 



(26) 



(27) 



(28) 



Note that Au is complex number, and that Qi and Q2 are complex conjugate with each other. On the other hand, 
Bi. is real as far as iliu is real, which is the case owing to our restriction to small background oscillations. 
In order to diagonalize the Hamiltonian ([27|) , we introduce the Bogoliubov transformation as 



with 



(29) 



Choosing 9k as 



af^ (t) = e'"^" cosh 0f^ , f3f^ (t) = e'"^^ sinh 0f^ . 

1^1 



tanh 29 f^ = 



Bi 



one can diagonalize the Hamiltonian into the form 



^k =^fc«fcW«fcW' 



(30) 



(31) 



(32) 



where use has been made of the Wronskian relation ( |26| ) to show that the coefficient of (|32[ ) is in fact given by Qfg . 
Therefore, the operator a^{t) defined in ( p9| ) does represent annihilation operator of physical quanta at time t. Using 
( pi| ) with (|8|) is is easy to show that 



cosh^0fe=^^^±^ 



2Vt 



sinh^ e 



k 



k 



^k -^k 
2f]i. 



(33) 



In the above treatment, we have rehed on the fact that amphtudes of the background oscillation is small so that fi^ 
is real. If we remove this assumption the system experiences instability in a certain period of time. It is not totally 
unstable system due to the nature of the instability and it would be interesting to work out such system. However, 
we do not engage such business in the present paper. Instead, we argue that it is important to take account of back 
reactions in really addressing such problem. It would even be the case that the quantum back reaction entirely cures 
the instability. 

Let us move on to the analysis of the eigenstates for the effective Hamiltonian. We define the vacua |0) and |0(t)) 
as 

«ofc|0>-0. afc(t)|0(i))=0, (34) 

and assume that the oscillator at i = is in the vacuum state |0), i.e., |0(t — 0)) ~ |0). It leads to the initial conditions 
for a^ (t) and /?^ (t) , 

l«fc(0)l = l, /3fc(0) = 0. (35) 

The physical quantum with momentum k is created (annihilated) by al{t) (afc(i)), because the Hamiltonian is 
diagonalized by them at each time t. As the system evolves, the state |0) is no longer the vacuum annihilated by 

«feW- . . . u^ . 

We rewrite the Bogoliubov transformation ( |29| ) into the form of unitary transformation: 

aofc=e-^=e-^fcWa^(t)e«fc«, 

(36) 
at^=e*^=e-^fe(*)at^(t)e°fc«, 

where 

^kit) = ^k (4°-fe " "fea^fc) • (37) 

From the definitions of vacua ( p4|) and of transformation (|3^), we have 

|0) =e"'^fc(*)|0(t)). (38) 

After a short calculation llioll , we obtain the vacuum state of the form 



|0) = TT , ^ exp 



—aUt)a^ . (t) 
2a* k^ ' -k^ ' 
k 



|0(t)). (39) 



The state is known as the squeezed state and widely used in quantum optics ||I^ . A significant feature of (|39| ) is that 
it has pairing correlation between back-to-back, namely fc and — fc, modes, and we will refer to it as the squeezed state 
of BCS type in this paper. In the next two sections, we will discuss, in detail, what are the characteristic features 
of particle production due to the parametric resonance mechanism. Earlier discussions of the squeezed state in the 
context of DCC may be found in [|l^,gO|. We also note that the particle production via the parametric resonance 
has been extensively discussed in the context of reheating in the infiationary universe ||2l|| - p3| , and a mechanism for 
producing superheavy dark matter |2J| . 

Before closing this section, we would like to comment on the relationship between our formalism and the one 
by Shtanov, Traschen, and Brandenberger p7| which can provide description of quantum systems with parametric 
resonances. Since we have adopted their formalism in our previous paper, we will try to compare and clarify the 
relation between various quantities that appear in both formalism. 

In short, there is a simple relationship between two formalisms: they start from the time-dependent operator (our 
af^{t)) which diagonalizes the Hamiltonian and then Bogoliubov transform it to the time- independent operator clqJ^- 
Therefore, such a correspondence is likely to exists. 

In order to make the correspondence more precise, we write down the equations obeyed by the coefficients of the 
Bogoliubov transformation in both formalisms. In STB formalism, they read 



/,\ ^^k n 2i r' O; dt' 

(40) 



2^fe 



o tj.\ ^k -2i f n,dt' 



On the other hand, the equations of aiu and /3i. in our formahsm are given by 



"fc(^) = ^^k 



2n 



(41) 



Above two sets of equations are consistent with each other provided that the phases of the coefficients satisfy the 
transformation rule: 

\-^Q^ + ^kQ^J J^ (42) 

Namely, the phase parts, which are undetermined in our present formalism, are chosen in a particular way in the STB 
formalism. Notice that the right-hand-sides of (^) are pure phases and hence they are consistent. 

IV. SINGLE PARTICLE DISTRIBUTIONS 

Let us discuss the qualitative features of the single particle distributions of particles produced to explore the 
parametric resonance mechanism. The multi pion and sigma state in our mechanism is given by 

\or®\or^\^), (43) 

where |0) stands for the BCS type squeezed state (p9|). The factorizcd form in (E3) stems from the linear approxima- 
tions described in Sec. II. 

The average number of quanta (n) u (t) with momentum k at time t is given by 

(n)fe(i) = (0|aj^(i)afe(i)|0). (44) 

We suppress in this section the isospin index, since the single particle distributions is independent of isospin. Given 
the squeezed state (p9|), it can be readily shown to be 

{n)kit)-W-^^^^ (45) 

K 

by using ( p3| ) with Bj^ defined in (Eq). 

Before presenting results of numerical computation, we discuss what would be the characteristic feature of the 
single particle distributions. It is natural to expect that an enhancement occurs due to the parametric resonance 
mechanism. With physical values of the parameters, we take the q parameters defined in (|^ and ( |l8| ) are given by 

q^ = 2.0 (^) , q, = 6.0 (^) . (46) 

Numerically, q-^ and q^ are less than ~ 0.1 and ~ 0.3, respectively, under the present approximation (^) < (^is^y. 
Therefore, we are working with narrow resonance approximation. One should keep in mind, however, that the narrow 
resonance approximation may not hold in the real world. 

At such small q parameter, the resonance occurs in the narrow bands at around the discrete values oi A, A — 
n^, (n = 1, 2, 3, ...). The first resonance takes place for pion at A^^ — 1 and it correspond to fc = 276 MeV. For sigma, 
the first resonance takes place at Ao- = 4 at zero momentum. The second resonance would be located at A^^ = 4 and 
at Act = 9 which imply k = 603 MeV and k = 692 MeV for pion and sigma, respectively. Since the parameter A is 
independent of ( — ), we expect that the peak position is stable against varying xo, as will be demonstrated below. 

We solve the Mathieu equation ( p3| ) subject to the boundary condition (p5|), and then compute {n)y,{t). In Fig. 
1-3 we plot the single particle momentum distributions of pion as a function of momentum k and dimensionless time 



z with background oscillation parameters — = 0.025, 0.05, and 0.1. We observe that a prominent peak exists at the 
right position for the first resonance in every three figures indicating the parametric resonance enhancement. In Fig. 
3 the momentum range |fc| < ^^O.Itti^ — rn\ in which an exponential instability exists (due to imaginary frequency) 
is cut off. 

In Fig.0, we give the single particle momentum distributions of sigma at ^ = 0.05. The first resonance, that is 
expected at fc = 0, is barely seen in the figure. The second resonance is invisible both in pion and sigma distributions. 

V. TWO PION CORRELATIONS 

We now discuss two particle distributions and correlations between pions. The quantum pion and sigma state is 
summarized in (Bsf) under the present approximations. We focus on the pion distributions of various charge states. 
Because of the factorized form of ( p9|) and (43) n-particle distributions of a can readily be obtained from that of tt" 
by an appropriate change of the parameters. It is true that the interpretation of our results on sigma particles in real 
experiments is not very transparent. Furthermore, when the a — *■ tttt coupling is turned on the a particles do affect 
the two pion correlations. Discussion of this point is, however, beyond the scope of this work. 

Two pion distributions are defined as 



^'^fc-fc'^ 



(V'l«2(i)afc(04(i)«i,WIV'), 



(47) 



Because of the factorized form of (^), there are no correlations between pion and sigma. Also, because of the 
factorized form of the squeezed state (p9h in momentum space, the nontrivial two particle correlations exist only 
between modes of identical momentum (h and fcV or between back-to-back momentum (k and — fc) configurations. 
One of the most important feature of the state (H3) is that it is the iso-singlet state. It comes from the fact that the 
frequency f2T is isospin singlet. As we discuss below, it will give us nontrivial constraints on two pion correlations of 
various charge combinations. At the same time it also generates the large charge- neutral fluctuations as discussed in 
Sec. I. 

In the computation of these quantities, we further recognize that the two-pion distributions at zero-momentum can 
not be obtained by taking the smooth limit fc — > of the expressions of either identical or back-to-back momentum 
configurations. It is because at_Q and al do not commute, whereas a^ and a , do commute for fe 7^ 0. Therefore, 
we have to compute three types of the two-pion distributions separately. 

We briefly describe some details of the computation of two pion distributions. It is convenient to use adjoint 
representation operator a\ for this purpose. The two pion distribution of the same momentum fc goes as follows: 



,t 



(^fc^fc) = («fc«fc«fc«fc 



= («fc«fc«fc«fc) + (a^afc) 



"fc 



"fc 



"fc 



(a_fca_fc«lfc«lfc> + |/5fc| 



((aLfca-fca-fc«Lfe) + 2(a_fca^^) j + l/^fel^ 
((«Lfc«-fc«Lfc«-fc) + 3(aLfca_A;) + 2 j + |/3fc| 



(48) 



"fc 



>-fe'^-fe>+3|/5fc| +2 +1/3;^ I , 



where the third equality stems from the relation 



^fc .t 



«fcl'^) = ;^-«-fel^)- 



(49) 



Similarly, (7r_fc7r_fc) is given as 



(^-k'^-k) = 



Pk 



4 



^fc^fc>+3|/3fcr + 2) + |/3fcr- (50) 



By solving the coupled equations (Hq) and (pO[), we obtain two pion distributions of j-th isospin component as 

Via analogous ways, one can compute the following expressions of the elements of nontrivial two pion correlations. 
For the same k ^ Q: 

Hk^'rk) = \f^k\'' ' 

^''Ik^'^k^lk^jk) = \f^k\'^ (for ' ^ j)' 

(52) 



"IkV^k^.k)^^ (forz^j) 



For opposite fc 7^ 0: 



(«Ifc"^-fc) = 0' 



(«Ife"'fc"!-fc^--fc) = l^fcl"* (for * 7^ j) 



(53) 

(«lfc^-fc«Uvfe) = l/^fcT (l/^fcT + 1 

By using these formulae, it is straightforward to compute the two pion distributions of definite isospin states. 
We transform them into the ones written by pion charge states, which are more suitable to gain insights for the 
experiments. We will describe in Appendix B the relationship between the adjoint representation pion operator (with 
which we have worked) and the charge-state operator, the (slightly modified) Horn-Silver representation 

k ^\ ^ ^ 

°r-;^(4 + »l). («) 

where a, and a, stand for the annihilation operators of tt* and tt" with momentum fc, respectively. We express the 
two pion distributions of various charge combinations in the form of the ratio R of them to the single pion distributions 
as defined by 

where p and q stand for ± and 0. They read: 
(a) identical momentum distribution: 

%,fe(+'+) = ^fc,fc(0'0)-2+-^, 

^ >^ (56) 

%,fc(+'") = 1- 
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(b) back-to-back momentum distribution: 

1 



Wfe (57) 



(c) zero-momentum distribution: 



%=o.fc=o(+' +) - %=o.fc=o(+' -) - 2 



(60) 



(61) 



(58) 
2 
^k=o.k=o(^^ 0) = 3 + — - — . 

As an another measure for two pion correlations, we define the correlation function in the following way: 

Then, it is straightforward to obtain the following results. 

(a) identical momentum correlations: 

C(^+,7r+) = C(4,4) = (n)|, 

(b) back-to-back momentum correlations: 

^H^^Zk^ = c(^l^-k^ - i^k) {i^^)k + ^) ^ 

(c) zero-momentum correlations: 

^(^Lo' '^fc^o) = <">fc=o (('^)fc=o + 1) ' (62) 

^(4=o''^Lo) = <")fc=o(2Mfc=o + l)- 

As will be discussed in Appendix B, the features of the correlation functions in (c) can be understood partly as a 
consequence of the isospin singlet nature of ( p3| ) . 

The identical and back-to-back momentum correlations are depicted in figures ra and o. In these figures, the 
amplitude parameter is taken as Xo/v — 0-05 and it is integrated over dimensionless time z from to 10, which 
corresponds to the period of time from i = to 6.5 fm. 

VI. POWER SPECTRUM AND POSSIBLE MECHANISMS FOR AMPLIFICATION OF 

LOW-MOMENTUM MODES 

Now we would like to discuss in some details the problem of the origin of the long lasting amplification of low- 
momentum pion modes seen in linear sigma model simulations, in particular in Ref. y. 

Let us first try to make a rough estimate of ^ in their simulation. It is shown in Fig. 2 of Q that (0^) is given 
by a constant superposed by an oscillatory component during the time scale of about 10-40 in units of a when large 
amplification takes place. Here a is the lattice constant and it is taken as a = (200MeV)^^ in W. Since the frequency 
of the oscillatory component reflect the pion mass not the sigma mass it cannot be used to estimate the amplitude 
of oscillation in sigma direction. We thus interpret (a quarter of) the constant component as a time and spatially 
averaged a oscillation. Since the low momentum modes are dominant at later times it may be utilized to estimate 
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the amplitude of uniform background a oscillations in our scheme. Using the constant component of ^ O.Ola ^ in 
((/)^) we obtain — ~ 0.16. It is a small number and indicates that our approximation is not totally absurd. But it is 
larger than (HhL)2 ^ 0.05, which implies that a new type of instability mentioned in Sec. Ill must also be excited. 

Since we confine ourselves into (^) < (^)^ in this paper, it is not possible to directly compare our parametric 
resonance mechanism to the Rajagopal-Wilczek simulation. Nevertheless, we attempt to make a bold comparison 
between these two by estimating the power spectrum, the corresponding quantity to the power calculated by them. 

Rajagopal-Wilczek used the definition of power spectrum pa 



PRw{k,t) = j^ 



^ e^''^^q^in,t) 



(63) 



where N indicates (one-dimensional) box size and they used A^ = 64 and (t>{n, t) denotes a single component of pion or 
a fields at lattice cite n. It has the dimension of [mass)'^ and hence the ordinates in their Fig.l must be understood 
to be plotted in units of a^^. 

We define the corresponding quantity, the power spectrum, in our framework. The natural definition is: 

Pours(fe,t)-(|Qfc(t)P>. (64) 

The expectation value (•) is to be evaluated by using the normalized squeezed state (p9). When the fields are confined 
into a finite volume V , Fourier integral in ( pl[ ) is modified to the discrete summation of the form 

e(a.,t)^^^e-fc^gfc(t), (65) 



y , 



Qj^it)^^ f d?xe^^^i{x,t). (66) 



or inversely, 

OiJt) = - 

Then, using V — (Na)^, our power spectrum (p4) is related with theirs via the following relation: 

PRw{k,t) = Pours{k,t)/a\ (67) 

It should be noticed that the number presented in W is PRwik, t)a? . Therefore, we must plot Poursik, t)/a to compare 
our power with theirs. 

We present our power spectra of pions and sigmas in Fig. 7-10. We have adjusted some parameters so that they 
accord with the values adopted in 0; wo = 87.4MeV, v = 92.5MeV, m^ = 135MeV, though its effect is tiny and does 
not affect our conclusions. We employ the amplitude parameter of background oscillation as Xq/v — 0.05 in order to 
avoid the instability. 

In Fig. 7 plotted is the power spectrum of pions at t = 0. One can see that there is a lot of power at relatively low 
momentum modes. One may feel curious of the fact that the power is not small at i = at which we set the boundary 
condition of no particle production. One can easily verify that it is natural because the latter condition implies that 

/O 1 1 



-4 , ^ , and Q '- 

^ V2S1 Vm,r ^ 

In Fig. 8 we present the time evolution of power of pion field at A: = 40MeV, which is the longest wavelength case 
depicted in Fig.l in B and it is where the largest amplification occurs. We observe that the average power is of order 
unity as we expect from the above estimation. It is a reasonable result because it is off resonance. 

With these parameters, the first resonance appears in fc = 268MeV and the amplification of the power in this 
momentum bin is depicted in Fig. 9. The Pours/o- is growing by a factor of 50 and reaches to ~ 20 which is similar 
order of magnitude with Rajagopal-Wilczek's. Of course, two numbers cannot be compared directly because of the 
difference involved between two computations including the initial conditions. Nevertheless, the qualitative feature of 
the growing power to the same order of magnitude may be an indication that the parametric resonance mechanism 
plays an important role in the long lasting amplification of low momentum modes. One may argue that our results 
and Rajagopal-Wilczek's are qualitatively different because the amplification occurs in low momentum modes in their 
simulation, whereas it happens at resonance in our case. However, if we have not ignored the nonlinearity it might 
have mediated the enhancement at resonance to other modes. 

In Fig. 10 plotted is the time evolution of power of sigma field at fc = 40MeV as in Fig. 8. The power of sigmas 
is an order of magnitude smaller than that of pions, in agreement with the feature obtained in Rajagopal-Wilczek's 
simulation. 
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We now summarize our present understanding on possible origin of the long lasting amplification of low-momentum 
pion modes seen in the linear sigma model simulations. We have suggested that the parametric resonance might be the 
cause. The semi quantitative features of the power spectrum including its order of magnitude and the rate of growth 
is not inconsistent with our proposal. We also pointed out that a new type of instability exists in our approximation 
scheme which, if real, would enhance the amplification of low-momentum modes. We, however, failed to achieve an 
intuitive understanding of the physical origin of the new instability. 

VII. CONCLUSION AND DISCUSSION 

We have discussed the pion production via the parametric resonance mechanism within the linear sigma model. 
In particular, stimulated by the feature of the numerical simulations of the linear sigma model, we focused on the 
scenario in which classical background oscillations of the sigma model fields are in the u direction. Assuming small 
amplitude oscillation which may be natural in the late stage of evolution of DCC, we have shown that one can ignore 
effects due to nonlinearity and quantum back reactions for sufficiently long time. Thanks to this fact, we were able 
to construct an explicit quantum pion (and sigma) state which allows us to calculate the two particle correlations as 
well as the single particle distributions. 

Our treatment, though under very restrictive assumptions, may be good enough to illuminate characteristic features 
of the parametric resonance mechanism. We formulated the quantum theory of the system with Mathieu instability 
on a more conventional basis of quantum field theory and clarified the relationship between our and the formalism 
given by Shtanov, Traschen and Brandenberger |]l7| . 

We have analyzed the computed single particle distributions and clarified the structure of narrow resonances charac- 
teristic to the parametric resonance mechanism. We then discussed the two pion correlations as a possible experimental 
probe for disoriented chiral condensate. Since the two pion correlations have unique characteristics, the back-to-back 
(in momentum space) correlations, it must give a clear signature which should merit the experimental hunting of 
DCC. In particular, it cannot be masked by the identical particle interference, the Hanbury-Brown-Twiss effect Pqj . 

Why two particle correlation? It is certainly far more difficult to measure compared with the multiplicity distri- 
butions, on which all the recent experimental search for DCC rely |g^,^. The global analysis using multiplicity 
distributions is powerful if large fractions of events are accompanied by the DCC domain formation. On the other 
hand, a different strategy is required for hunting if DCC is a rare phenomenon. 

We have discussed possible origin of the long lasting amplification of low-momentum pion modes in the last section. 
Our discussion cannot be a complete one, but we hope that it stimulates the readers' interests in this problem. 

We should mention about the limitation and the drawback in our treatment in this paper. We ignored the quantum 
back reaction and the instability of a meson which should exist in the real world. It is the perfectly legitimate 
approximation if the amplitudes of background a oscillations are really small. However, it is possible that the 
amplitudes are sometimes large because of the prevailing thermal fiuctuations in the initial stage. Since we are dealing 
with the system in which the coupling is really strong, most probably, the peak in the single particle distributions 
will go away after the quantum back reaction is taken into account. 

Then, the key question is that if anything in our results remain valid after the quantum back reaction is taken 
into account. We argue that the answer is yes; it is the characteristic features of the two pion correlations which are 
discussed in detail in Sec. V. We now engage in a computation to verify our expectation. 

How about the a — *■ tttt coupling? It is clear that it also tends to smear out the resonance peaks and, more 
importantly, may obscure the signature of the back-to-back correlations. Again we need more elaborate treatment 
which includes the instability of a to make definitive statement about how much the signature survive in the case with 
a — > TTTT coupling. The formalisms which may be suitable for such analysis have been investigated in detail p9| , p0|] . 

There are also some recent proposals p^p^ that the hadronic medium effects may induce the similar back-to-back 
correlations in momentum space, which would mimic the signature of DCC discussed in this paper. One would 
hope that it should be possible to find observational features which discriminate these two mechanisms. The task is, 
however, beyond the scope of this paper. 

ACKNOWLEDGMENTS 

During the long term investigation we benefited from numerous conversations in various stages of this work with 
Masayuki Asakawa, Dan Boyanovsky, Berndt Miiller, Krishna Rajagopal, and Valery Rubakov. In particular we 
thank Krishna for his informative and kind correspondences in comparing power spectra obtained by our and their 
calculations. 



13 



The research of HM is partly supported by the Grant-in- Aid for International Scientific Research No. 09045036, 
Inter-University Cooperative Research, and the Grant-in- Aid for Scientific Research in Priority Areas No. 11127213, 
Ministry of Education, Science, Sports and Culture of Japan. 

APPENDIX A 

We discuss in this Appendix some aspects of multiparticle correlations among particles involved in a normalized 
single-mode squeezed state 

IWA = i ^-^^+10), (A.l) 



v/(0|e^-^-e^^+|0) 

-(a)^ and - 
N = a}a, form a closed algebra. 



where Kj^ = 77 1*^) ^-nd K^ = :^(o.) , and 7 denotes a complex number. The operators Kj^ and K^^ together with 



[if_,if+]=iV+i, [iV,K±]=±2K±, (A.2) 

which will play a role in our following computation. 

For convenience in the systematic treatment, we employ the generating function formalism developed by Koba, 
Nielson, and Olesen, |^^ and by Koba [Q. The generating function is defined by 

C30 

F\h\^Y.^\ + hTP,,, (A.3) 

using the multiplicity distribution P„ = |(7i|?/;)p where 

\n) = i^ |0). (A.4) 

\/n\ 

The generating function can generate whole set of various moments when expanded in various different manner: 

n=0 

(A.5) 

00 , /; 



exp 



h' 



.k=0 



In these equations, F{k) denotes 

F^^) = {n{n -l){n-2)---{n~k + 1)), (A.6) 

where (• • •) implies the average over by the multiplicity distribution P„; (O) = X]n=o OnPn- As is familiar in cluster 
expansion in statistical mechanics, C*-*^' represents the correlations 

C(i) = (n), 

C(2) = (n(n-l)>-(n)2, (A.7) 

C(3) = (n(n - l)(n - 2)) - 3(n(n - l))(n) - (n)3, 

and so on. 

The expression of the generating function in a form of operator expectation value is give by Koba p4| ; 

F[h] = (^1 : e^'^'^ : |^), (A.8) 
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where :: indicates to take the normal ordering. It proves to be a very useful formula for our purpose. Toward 
calculating F[h], we define the quantity Am as Am = {4>\{a'')"^{a)"^\'ip). Using the algebra applied to the vacuum 
state 



m „7-R'-t- 



la"%e''^+J |0) =7a™^iaTe'^-f^+|0), 



(A.9) 



one can derive the recursion relation among A„ 

Ao =1, 
Ai = (n), 



(A.IO) 



Am+2 = (n) 



(2m + 3)Am+i + (to + 1)2A„ 



(to > 0). 



One can show, by direct computation, that 



We can convert the recursion relation ( A.IOJ) into the differential equation obeyed by the generating function 



(A.ll) 



1 . d^ 



d 



{h' +2h--^)--r+i{h+ 1)— + 1 



dh 



dh 



F[h] = 0. 



We make a change of variable 



i = 



[l + hf 



l + (n 
to transform the differential equation into the standard form 

^ d- ,1 ^ d 1 

The solution to the equation subject to the boundary conditions 



F[t] = 0. 



(A.12) 



(A.13) 



(A.14) 



is uniquely given by 



F[t] 



nt] \,-_ 



F[t] 



!+("> 



-1, 



1 + (n) 
2 ' 



(A.15) 



V(l + (n»(l-t) 



(A.16) 



The boundary conditions correspond to Aq — 1, Ai ~ (n), and A2 = {n){3{n) + 1). Thus, we have obtained the 
generating function as 



F[h] 



The correlation moments can be easily computed as follows: 



(A.17) 
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C(2) =(n)(2(n) + l), 
C(3) =2(n)2(4(n) + l), 



(A.18) 



Since the generating function (A. 17) is the function of (1 + ft,)^, the muhiphcity distribution P„ for odd n vanishes. 
That for even n is given by 



P2k = 



(2fc)! 



1 



22'=(fc!)2yrT7;iy Vl + (n) 



(A.19) 



which is n othing but the negative binomial distribution famihar in hadronic multiparticle phenomenology. The 
expression ( A.19 ) reproduces the results obtained by Yoshimura |23| who also calculated all the off-diagonal elements 
of the density matrix. 



For completeness, we calculate the KNO scaling function for the multiplicity distribution ( A.19 ). It is obvious, 
from the expression of the correlation moments, that the squeezed state has "long range" correlations 



(jk 



finite. 



(A.20) 



the multiplicity distribution obeys the KNO scaling behavior. Either by solving the moment problem or by using the 
explicit form (A.19), one can show that 



^{z) 



lim(„)^oo('^)-Pn 



Ae--/2 



(A.21) 



APPENDIX B 



We derive relations among two-pion correlation functions of various different charge states which follow from isospin 
invariance. We restrict ourselves into the case with iso-singlet state which is of concern to us. We also confine ourselves 
into the case of zero momentum pions; the treatment for nonzero momentum pions is different and is much more 
involved. 

We work with the (slightly modified) Horn-Silver pq] representation of the isospin generator 



T — a^ Tij a-' , 



(i,j = +,-,0) 



where 



Tl 



1 

71 



[0 





1' 








1 


1 


1 






T2 



1 

72 









—i 








i 


i 


—i 






T2 



1 





0] 





-1 















The explicit form of T is as follows: 

V2 
V2 



T^3 __ 

V2 



i(-)'a(o) 






a(+)^a(+)-a(-)^a(-) 



(B.l) 



(B.2) 



(B.3) 
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where a^"*"-*, a^°-*, and a^ ^ are the annihilation operators of 7r+,7r°, and tt , respectively, as defined in (pi]). The 
relation between the charge eigenstate operators and the adjoint representation operators are given in (|54|). 

Let \ip) be the iso-singlet state. Then, it is annihilated by the isospin operator, T\tjj) — 0. From this, it follows that 

(iPle'^-^a'^^a^li;) = (i^|a'=^a'|i^), (B.4) 

for any k, I. This is the relation that should be obeyed for all orders in ^. If we take the first-order term in ^, we 
obtain the sum rule 

(V'lfa'^^a'lV^) =0 {for any k, I). (B.5) 

Due to isospin invariance, it is easy to show that (7t+) = (tt^) = (tt*^). Then, the T^ (or T^) sum rule produces the 
relations like 

(^+^0) = (tt-tt"), (B.6) 

which also follows from the T^ sum rule. The useful relations comes from T^-sum rule 

(7r+7r+> = (tt+tt-) = {-r-t^-). (B.7) 

In these equations, {n^n^) implies the two-pion distribution of charge state a and /3. 
The two-pion correlation function may be defined as 

C(7r+,7r+) = (7r+7r+) - (7r+) - (7r+)2, 

C{7T+,Tr°) ^ (tt+tt") ~(Tr+)(Tr°), (B.8) 

etc. 

Then, we obtain the isospin sum rule, 

C{TT+,TT-)^CiTT+,TT+) + {TT), (B.9) 

where (tt) = {tt+) = (tt^) = (tt*^). Since (tt) is positive, the inequality among two-pion correlations follows; 

C{tt+,tt-)>C{tt+,tt+) (B.IO) 

It may be surprising that the unlike-sign correlation is stronger than like-sign correlation. But it is the consequence 
of the isospin invariance with iso-singlet nature of the multipion state. 
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FIG. 1. The single pion momentum distribution as a function of the scaled time z — i(mo-i + 7r). The amplitude of oscillation 
is taken as Xo/v = 0.025. 

FIG. 2. The single pion momentum distribution as a function of the scaled time z = ^(mct + n). The amplitude of oscillation 
is taken as Xo/'v = 0.05. 

FIG. 3. The single pion momentum distribution as a function of the scaled time z = i(mo-i-|-7r). The amplitude of oscillation 
is taken as Xo/^ = 0.1. 

FIG. 4. The single sigma momentum distribution as a function of the scaled time z = ^{nicrt + n). The amplitude of 
oscillation is taken as Xo/" ~ 0.05. 
FIG. 5. The pion identical momentum correlation function is depicted. The amplitude of oscillation is taken as xo/'" = 0.05. 

FIG. 6. The pion back-to-back momentum correlation function is depicted. The amplitude of oscillation is taken as 

Xo/v = 0.05. 

FIG. 7. The power spectrum of pion at t = is depicted. Due to the boundary condition of no particle production at i = 0, 
the power spectrum is not small. 

FIG. 8. The time evolution of power spectrum of pion field at fc = 40MeV is presented. It corresponds to the longest 
wavelength bin in Fig. 1 of Rajagopal and Wilczek's paper. 

FIG. 9. The time evolution of power spectrum of pion field at momentum k — 268MeV, which corresponds to the first 
resonance band, is drawn. 

FIG. 10. The time evolution of power spectrum of sigma field at fc = 40MeV is depicted. 
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